Entropy production and time asymmetry in nonequilibrium fluctuations 
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The time-reversal symmetry of nonequilibrium fluctuations is experimentally investigated in two 
out-of-equilibrium systems namely, a Brownian particle in a trap moving at constant speed and an 
electric circuit with an imposed mean current. The dynamical randomness of their nonequilibrium 
fluctuations is characterized in terms of the standard and time-reversed entropies per unit time of 
dynamical systems theory. We present experimental results showing that their difference equals the 
thermodynamic entropy production in units of Boltzmann's constant. 
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Newton's equations ruling the motion of particles in 
matter are known to be time-reversal symmetric. Yet, 
macroscopic processes present irreversible behavior in 
which entropy is produced according to the second law of 
thermodynamics. Recent works suggest that this thermo- 
dynamic time asymmetry could be understood in terms 
similar as those used for other symmetry breaking phe- 
nomena in condensed matter physics. The breaking of 
time-reversal symmetry should concern the fluctuations 
in systems driven out of equilibrium. These fluctuations 
may be described in terms of the probabilities weighting 
the different possible trajectories of the systems. Albeit 
the time-reversal symmetry of the microscopic Newto- 
nian dynamics says that each trajectory corresponds to 
a time-reversed one, it turns out that distinct forward 
and backward trajectories may have different probability 
weights if the system is out of equilibrium. For example, 
the probability for a driven Brownian particle of having 
a trajectory from a point A to a point B is different of 
having the same reverse trajectory from B to A. 

This important observation can be further elaborated 
to establish a connection with the entropy production. 
We consider the paths or histories z = (zq, z\, Z2, Zn-i) 
obtained by sampling the trajectories z(t) at regular time 
intervals r. The probability weight of a typical path is 
known to decay as 

P + (z , Zi, Z 2 , ... ,Zn-i) ~ exp(-7lT/l) (1) 

as the number n of time intervals increases [3, 0, [H, IH . 
The decay rate h is called the entropy per unit time and it 
characterizes the temporal disorder, i.e., dynamical ran- 
domness, in both deterministic dynamical systems and 
stochastic processes [H, H, [H, Hj]. We can compare (HJ) 
with the probability weight of the time-reversed path 
z^ = (z„_i, Z2, zi, Zq) in the nonequilibrium system 
with reversed driving constraints (denoted by the minus 
sign): 

P-(z n -i, ...,Z2,Zi,Zq) ~ exp(-nr/i R ) . (2) 



It can be shown that, out of equilibrium, the probabilities 
of the time-reversed paths decay faster than the proba- 
bilities of the paths themselves [Bj]. We may interpret 
this as a breaking of the time-reversal symmetry in the 
invariant probability distribution describing the nonequi- 
librium steady state, the fundamental underlying Newto- 
nian dynamics still being time-reversal symmetric. The 
decay rate h K in Eq. ((2|) is called the time-reversed en- 
tropy per unit time and characterizes the dynamical ran- 
domness of the time-reversed paths [1, In the case 
of Markovian stochastic processes with discrete fluctuat- 
ing variables, the difference between both quantities /i R 
and h gives the entropy production of irreversible ther- 
modynamics 0, @> S Hi- A closely related result has 
been obtained for the work dissipated in transient time- 
dependent systems [||. However, many experimental sys- 
tems have continuous fluctuating variables and evolve in 
nonequilibrium steady states. Therefore, we may wonder 
how to measure dynamical randomness in such systems 
of experimental interest and whether the time asymme- 
try of this property can be experimentally detected and 
related to the thermodynamic entropy production. 

In the present Letter, we provide experimental ev- 
idence for the aforementioned time asymmetry in two 
nonequilibrium systems, namely, a driven Brownian mo- 
tion and a fluctuating electric circuit. For this pur- 
pose, the decay rates h and h n are considered as so- 
called (e, r)-entropies per unit time, which characterize 
dynamical randomness in continuous-variable stochastic 
processes These entropies per unit time can be ob- 
tained by applying to the present stochastic systems a 
method originally proposed for the study of deterministic 
dynamical systems [H, El [H . Thanks to this method, the 
(e, r)-entropies per unit time of the paths and the corre- 
sponding time-reversed paths can be evaluated from two 
long time series measured with sufficient temporal and 
spatial resolutions, in two similar runs but one driven 
with an opposite nonequilibrium constraint. The experi- 
ment thus consists in recording a pair of long time series 
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in each system. The dissipated heat and thermodynamic 
entropy production are thus given by the difference be- 
tween the two (e, r)-entropies per unit time. 

The first system is a Brownian particle dragged by an 
optical tweezer, which is composed by a large aperture 
microscope objective (x63, 1.3) and by an infrared laser 
beam with a wavelength of 980 nm and a power of 20 
mW on the focal plane. The trapped polystyrene parti- 
cle has a diameter of 2 fxm and is suspended in a 20% 
glycerol- water solution. The particle is trapped at 20 /im 
from the bottom plate of the cell which is 200 /jm thick. 
The detection of the particle position x t is done using 
a He-Ne laser and an interferometric technique (To| . In 
order to apply a shear to the trapped particle, the cell 
is moved with a feedback-controlled piezo which insures 
a perfect linearity of displacement. The motion of the 
dragged particle is overdamped and can be modeled as 
the Langevin equation 



dx t 

a ^r = F{xt 



ut) + , 



(3) 



where a is the viscous friction coefficient, F = — d x V is 
the force exerted by the potential V — kx 2 /2 of the laser 
trap moving at constant velocity u, and £(t) a Gaussian 
white noise The stiffness of the potential is k = 

9.62 10 -6 kg s -2 . The relaxation time is = a/k = 
3.05 10" 3 s. 

The second system is an electric circuit driven out of 
equilibrium by a current source which imposes the mean 
current / 11211. The current fluctuates in the circuit be- 



cause of the intrinsic Nyquist thermal noise [ll|. The 
electric circuit is composed of a capacitance C — 278 pF 
in parallel with a resistance R — 9.22 Mil so that the 
time constant of the circuit is tj? = RC — 2.56 10~ 3 s. 
This electric circuit and the dragged Brownian particle, 
although physically different, are known to be formally 
equivalent by the correspondence a <-> R, k <-> 1/C and 
u <-> I while the particle position x t corresponds to the 



charge qt inside the resistor at time t [ll|, LL2J . The vari- 
ables Xt and qt are acquired at a sampling frequency 
1/t = 8192 Hz. 

In both experiments, the temperature is T = 298 K. 

In order to fix the ideas, we describe our method in 
the case of the dragged Brownian particle. The heat dis- 
sipated along a random trajectory during a time interval 
t is given by 11. 13j 
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dx t i 
~dt 



T F(x v - ut 1 ) dt' . 



(4) 



After a long enough time, the system reaches a nonequi- 
librium steady state, in which the entropy production is 
related to the mean value of the dissipated heat according 
to 



a\S = 1 d(Q t ) 
dt T dt 



Our aim is to show that one can extract the heat dis- 
sipated along a fluctuating path by comparing the prob- 
ability of this path, with the one of the time-reversed 
path having also reversed the displacement of the po- 
tential, i.e., u — ► — u. We first make the change to the 
frame comoving with the minimum of the potential so 
that z = x — ut. After initial transients, the system will 
reach a steady state characterized by a stationary prob- 
ability distribution. As we are interested in the proba- 
bility of a given succession of states corresponding to a 
discretization of the signal at small time intervals r, a 
multi-time random variable is defined according to Z = 
[Z(to),Z(to + r), . . . , Z(t + nr — t)] which corresponds 
to the signal during the time period t — to — nr. For 
a stationary process their distribution do not depend on 
the initial time t$. From the point of view of probability 
theory, the process is defined by the rt-time joint proba- 
bilities P a (z; dz, t, n) = Pr{z < Z < z+dz; a} = p a (z)dz, 
where p a (z) is the probability density for Z to take the 
value z = (zo, Z\, . . . , at times to + ir for a nonequi- 

librium driving a = u/\u\ = ±1. Since the process is 
Markovian, the joint probabilities can be decomposed 
into the products of the Green functions G{zi, z^_i; r)dzi 
for i = 1, ... ,n. G{z, zo; t) gives the probability density 
for the position to be z at time t given that the initial po- 
sition was zq [3, EH- To extract the dissipation occurring 
along a single trajectory, one has to look at the ratio of 
the probability of the forward path over the probability 
of the reversed path having also reversed the displace- 
ment of the potential. Indeed, taking the logarithm of 
this ratio and the continuous limit r — > 0, n — > oo with 
tit = t, we find 
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P+(z; dz, r, n) 



I3u / F(zf) dt'-p V(z t )-V(z Q ) 



au 



(5) 



P_(z R ;dz, r,n) 

(6) 

which is exactly the heat Qt in Eq. ([4]) expressed in 
the z variable and multiplied by the inverse temperature 
(3 = (/cbT) -1 . We notice that, alone, the first term gives 
the work exerted by the trap 

nana. 

Relations similar to Eq. ([6]) have been obtained for the 
distribution of the work done on a time-dependent sys- 
tem 17 , HI and for Boltzmann's entropy production [3] • 
We emphasize that Eq. ([6]) also holds for anharmonic po- 
tentials V and that the reversal of u is essential to get 
the dissipated heat from the way the path probabilities 
P + and P_ differ. 

Now, due to the continuous nature in time and in 
space of the process, one has to consider (e, r) quantities, 
i.e. quantities defined on cells of size e and measured at 
time intervals r. Therefore, we introduce the probability 
P + (Z m ;e,r, n) for the path to remain within a distance 
e of some reference path Z m , made of n successive posi- 
tions of the Brownian particle observed at time intervals 
r for the forward process. The probability is obtained by 
searching for the recurrences of M such reference paths 
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FIG. 1: Time series of typical paths z(t) for the Brownian 
particle in the optical trap moving at the velocity u for the 
forward process (upper curve) and — u for the reversed process 
(lower curve) with u = 4.24 10~ 6 m/s. 



or patterns in the time series. Next, we also introduce 
the probability P_(Z^; e, r, n) for a reversed path of the 
reversed process to remain within a distance e of the refer- 
ence path Z m (of the forward process) during n successive 
positions. According to a numerical procedure proposed 
by Grassberger, Procaccia and others [2, Q the entropy 
per unit time can be estimated by the linear growth of 
the mean 'pattern entropy' defined as 



1 M 

H(e,T,n) = -— lnP+{Z m ;e,T 



By similarity, we introduce 



(7) 



1 M 



(8) 



for the reversed process. The (e,r)-entropies per unit 
time, h(e, r) and h R (e,r), are defined by the linear 
growth of these mean pattern entropies as a function of 
the time m [H, 0, 0] . In the nonequilibrium steady state, 
the thermodynamic entropy production should thus be 
given by the difference between these two quantities: 



h(e,r)] 



(9) 



It is important to note that the probabilities of the re- 
versed paths are averaged over the paths of the forward 
process in order for Eq. © to hold. The entropy pro- 
duction is thus expressed as the difference of two usually 
very large quantities which increase with the scaling law 
e~ 2 for e,r going to zero [3, 01. Nevertheless, their dif- 
ference remains finite and gives the entropy production in 
terms of the time asymmetry of the dynamical random- 
ness characterized by the (e,r)-entropies per unit time. 
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FIG. 2: (a) Entropy production of the Brownian particle 
versus the driving speed u. The solid line is given by Eq. 
@. (b) Entropy production of the RC electric circuit ver- 
sus the injected current /. The solid line is the Joule law, 
diS/dt = RI 2 /T. The dots are the results of Eq. ©. 



In order to test experimentally that entropy produc- 
tion is related to this time asymmetry according to Eq. 
([9]), we have analyzed for specific values of \u\ or |/|, a pair 
of time series up to 2 10 7 points each, one corresponding 
to the forward process and the other corresponding to 
the reversed process, having first discarded the transient 
evolution. Figure Q] depicts examples of paths z(t) for the 
Brownian particle in a moving optical trap. 

For different values of e between 5.6-11.2 nm [2lj |. the 
mean pattern entropy ([7]) is calculated with the distance 
defined by taking the maximum among the deviations 
\Z(t) — Z m (t)\ with respect to some reference path Z m 
for the times t = 0, r, . . . , (n — l)r. The forward entropy 
per unit time h(e, t) is evaluated from the linear growth 
of the mean pattern entropy ([7]) with the time nr. The 
backward entropy per unit time h R (e, t) is obtained sim- 
ilarly from the time-reversed pattern entropy ([8]). The 
difference of the two dynamical entropies is depicted as 
in Fig. (2^.. The good agreement with the entropy pro- 
duction is the experimental evidence that this latter 
is indeed related to the time asymmetry of dynamical 
randomness as predicted by Eq. ([5]). 
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On the other hand, we have analyzed by the same 
method the time series of the RC electric circuit. We see 
in Fig. [2h that the entropy production obtained from 
the time series analysis of the RC circuit agrees very well 
with the known Joule law, which is a further confirmation 
ofEq. ©. 

We also tested the possibility to extract the heat ((3]) 
dissipated along a single stochastic path by searching for 
the recurrences in the time series according to Eq. |6|). 
A randomly selected path as well as the corresponding 
heat dissipated are plotted in Fig. [3l We find a very 
good agreement so that the relation ^ is also verified for 
single paths. In this case, the heat exchanged between 
the particle and the surrounding fluid can be positive or 
negative because of the molecular fluctuations. It is only 
by averaging over the forward process that the dissipated 
heat takes the positive value depicted in Fig. O 

In conclusion, we measured the entropy production by 
searching the recurrences of trajectories in the fluctuat- 
ing dynamics of two nonequilibrium processes. The ex- 
periments we performed consisted in the recording of two 
long time series. The first one corresponds to a forward 
experiment while the other is measured from the same 
experimental setup except that the sign of the constraint 
driving the system out of equilibrium has been reversed. 
From these two time series, we are able to compute two 
dynamical entropies, the difference of which gives the en- 
tropy production. Moreover, we tested the possibility to 
extract the dissipated heat along a single random path. 
This shows that the entropy production arises from the 
breaking of the time-reversal symmetry in the probability 
distribution of the statistical description of the nonequi- 
librium steady state. Since the decay rates of the multi- 
time probabilities of the forward and reversed paths char- 
acterize their dynamical randomness, the present results 
show that the thermodynamic entropy production finds 
its origin in the time asymmetry of the dynamical ran- 
domness. 
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